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Abstract 

I show that P(2,3) has an embedding P(2,3) S T C P(4,5,6,9) 
whose image T is contained in a quasismooth K3 hypersurface A24 C 
P(4, 5, 6,9). The pair T C X24 unprojects to the codimension 4 K3 
surface Y C P(4, 5, 5, 6, 7, 8, 9) with 

Basket = 4), ±(2, 3), i(4, 5)]. 

Numerator = 1 - t 12 - t 13 - 2t 14 - 2t 15 - 2t 16 - t 11 

+ t 19 + 2t 20 + it 21 + At 2 + 3t 23 +■■■ 

(Altmok4(lll) in the Magma K3 database). The local coordinates at 
the third centre P3 = |(2, 3) of Y are of weight 7 and 8 (rather than 
2 and 3), so both are eliminated by the projection from P3. Together 
with other examples, this gives substance to Type IV unprojections. 
Several more cases of Type IV unprojections are known up to codi- 
mension 5 or 6. The paper also contains some Magma programming 
routines suitable as exercises for babies. 



1 Introduction 

Unprojection makes new Gorenstein rings out of old. The idea is to use 
the adjunction formula in Serre-Grothendieck duality to construct rational 
functions with poles along a given divisor. Adjoining these functions is 
an analog of Castelnuovo's contractibility criterion in minimal models of 



surfaces. Type I unprojections were introduced by Kustin and Miller [KM] 



and later by Papadakis and Reid [PR]. They also figure prominently in 



birational geometry, see Corti, Pukhlikov and Reid [CPR], and |Kj], and in 



that context can be traced back to work on the Cremona group in the late 



19th century (Geiser involutions, see CPR], 2.6.3.). Unprojections become 



progressively more complicated as the divisor to be unprojected gets further 
from projectively Gorenstein, so that several unprojection variables have to 



be adjoined at the same time. Type II projections first occurred in [CPR], 
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4.10-4.12 and 7.3 (in which context they are related to Bertini involutions), 



and are studied in more detail in [Ki], Section 9. The typical case of Type III 
projections is the classical projection of a Fano 3-fold V 2 g~ 2 C P 9+1 from a 
line (compare JK|, Section 9 and Example 9.16). 

This paper introduces Type IV unprojections. I don't intend to give a 
formal definition here. The defining feature of Type IV is that the coordinate 
ring of the divisor to be unprojected is a ring A = C[T] whose normalisation 
A is Gorenstein, and such that the quotient A/ A is an Artinian module 
needing 2 linearly independent generators over A. An ideal case would be 

A = C[r] := C[u 2 ,uv,v 2 ,u 3 ,u 2 v,uv 2 ,v 3 }, with A = C[u,v]. (1.1) 

Here A needs l,u,v as generators over A, and u,v are linearly indepen- 
dent. I do not treat this case here, because any example with it would have 
codimension > 4 unprojecting to codimension > 7; see however Section |6| 

The Type IV unprojections studied in this paper add 3 new variables 
as generators to the ambient ring, so potentially increase the minimal num- 
ber of generators by 3 (hence also the codimension). But this increase is 
masked if some of the old variables become expressible in terms of the new; 
this frequently happens in small codimension. The K3 database in Magma 



(export 2.8, [Ma]) reports 5 families of K3 surfaces in codimension 4 hav- 



ing a Type IV unprojection to codimension 1 (see 7^ and compare | ABR| ); 



there are also 6 families in codimension 5 having a Type IV unprojection to 
codimension 2 and several similar things in codimension 6. 



2 Embedding of T and quasismoothness of X24 

Write u, v for coordinates on P(2, 3) and x,y,z,t for P(4, 5,6,9). Up to 
change of coordinates, the general map i : P(2, 3) — > P(4, 5, 6, 9) is 

i : (u, v) i-> (u 2 ,uv, u 3 + v 2 , u 3 v + v 3 ). (2.1) 

Choosing the 4th entry t = u 3 v + v 3 (so that t = vz) tidies up later calcu- 
lations. It is an embedding because monomials in (u 2 , uv , u 3 + v 2 ,u 3 v + v 3 ) 
span the vector space H°(P(2, 3), 0(i)) for alH 3> (in fact for every i > 13, 
see [74]). 

The image F = i(P(2,3)) is defined by 

/i8 : x 3 z + xyt — z 3 + t 2 , (721 : x 3 t — xyz 2 + y 3 z, 

gig: x 2 yz - xzt + y 2 t, q 22 : -xt 2 +y 2 z 2 , (2.2) 

h 2 o : x 5 - x 2 z 2 + 2xy 2 z - y 4 , q 2 3 : -x 2 zt + yz 3 - yt 2 
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(see |7.1| ). For reasons of degree, only the first 3 of these can take part in the 
equation of X24: set 

F 2 4 = xh 2 o + ygi9 + zfi 8 . (2.3) 

One checks by brute force (see ^J) that F defines a quasismooth hyper- 
surface X 24 C P(4,5,6,9). 



3 Resolution of C[u, v] as a C[x,y,z,t] module 

Write A = C[x,y, z,t] for the homogeneous coordinate ring of P(4, 5,6,9) 
and C[r] = A/Ir for that of T C P(4, 5, 6, 9); by construction, its normalisa- 
tion is the coordinate ring C[u, v] of P(2, 3). Thus the embedding i of ( |2.l| ) 
makes C[it, u] a module over the polynomial ring C[x,y,z, t], with 

x ■ 1 = u 2 , v ■ 1 = uv, z ■ 1 — x ■ u = v 2 , x ■ u = n 3 , 

(3.1) 

x ■ v = u 2 v, y ■ v = uv 2 , z ■ v — x 2 ■ 1 = v 3 , etc. 

This already shows that C[u, v] is generated as a C[x, y, z, t]-module by the 
3 elements (l,u,v). The relations between these generators are 



where 



(y, u, 1)M = 0, (3.2) 



x —y —z —tO 
M= I y z -xz t x 2 ) (3.3) 

—x 2 t z 2 —yz —xz + y 2 



(see 7.2). The 3x6 matrix M is homogeneous with entries of weight 

'4 5 6 9 8 7' 

5 6 7 10 9 8 I (3.4) 



,7 8 9 12 11 10, 



and satisfies 



MJ*M = 0, where J = q 3 ) • ( 3 - 5 ) 

Now C[u, v] is a Gorenstein module of codim 2 over C[x, y, z, t], and one 
checks that its resolution is 

C[u, v]< — L ^- L x J^- L 2 *- 0, (3.6) 
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where 



L = A©A(-2)©A(-3), 

La = A(-7) © A(-8) A(-9) A(-12) © A(-ll) © A(-10), 
L 2 = A(-16) © A(-17) © A(-19). 

4 Maps between complexes 

I find the unprojection variables by comparing the resolutions of C[X] = 
A/{F2a) and C[u, v]. For this, consider first the diagram including part of 
the resolution of C[T]: 



C[X] < 


- A 


«- A(-24) 


I 




I 


c[r] < 


- A 




n 


n 


1 




- L 





(4.1) 







Here Jfi = A(— 18) © ■ ■ ■ © (—23) is the free module corresponding to the 
6 generators ( |2.2j ) of the ideal Jr- The first downarrow A (—24) — > JTi 
is the matrix (z, y, x, 0, 0, 0) that expresses J24 as the combination Q2.3] ) of 
/i8) 519) ^20- For any / G Jp, the column vector (0, 0, /) G ker{Jo — * u]}, 
and is therefore hit by some element of L\. The second downarrow is the 
6x6 matrix expressing the 6 generators ( |2.2| ) of i£i as linear combinations 
of the 6 columns of M. This matrix is presumably implicit in the Grobner 
basis calculation that eliminates u, v. Because J24 is a linear combination 



of only the first three generators (2.3), I only need the corresponding three 
columns: 



I = MN, (4.2) 



Jl8 319 h 



20/ 
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where 



N 



' 
— xz 
t + xy 
— z 










-x 2 
t 



x 2 y + xt\ 
—x 3 


-x 2 
-xy 

-y 2 



(4.3) 



Composing gives 



C[X] * 


- A <r 


- A(-24) <- 


- 




i 


n 


i 






C[u,v] <• 


- L <- 






- 



(4.4) 



where the downarrow takes the basis -F24 to 



/ x 3 y + x 2 t\ 

4 2 
—x — xz 

xyz + zt 

—x 3 — z 2 

„2„ 




—2xy 
\-xy 2 + ytj 



G Li(24). 



(4.5) 



This allows me to start writing down the unprojection Y of T in X. The 
theory is similar to that of Kustin and Miller [ KM[ , and Papadakis and Reid 
|PRH and [Ki]: the adjunction formula ojy = £xtx(Or,u)x) gives an exact 



sequence 



— > lux —>■ T~(-omo x {1r,ux) — >• — > 0. 



(4.6) 



Note that I abuse notation by writing X, T in place of the punctured affine 
cones over them. The modules and Horns between them are really the Serre 
modules of sheaves on X, F. For example, T = P 1 , but as P(2, 3) it has 
cur = C(— 5), since on the punctured affine cone over T, the dualising sheaf 
twisted by 0(5) is based by du A dv. 

Moreover since C[r] C[u, v] is an isomorphism outside the origin, the 
dualising module Wcrr] = w c[u,i)] is isomorphic to C[u,v], so is free of rank 1 
as a C[u, v] module; however, as a module over C[r] or C[X] or C[x,y, z,t], 
it needs 3 generators, corresponding to l,u, v. I used the fact that C[u, v] is 
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a Gorenstein module to derive the resolution ( |3,6| ); it is best not to enquire 
too closely about the resolution of C[T]. 

Now the homogeneous coordinate ring C[Y] of the unprojected vari- 
ety Y is obtained from C[X] by adjoining rational functions so,s±,S2 with 
poles along T, viewed as homomorphisms Zp — > = Ox- Here sq G 
7^om(Tr, is the rational form with pole along T whose Poincare 



residue is a basis of or (5) 



O 



P(2,3)- 



As in the proof of Castelnuovo's 



criterion, the graded ring with sq adjoined already contracts T to a point. 
The elements s\,S2 that map to u, v times the generator of ct>r(5) under 
TComo x (Ir,^x) -» <^r are required to make the image projectively normal. 

I choose so to map to 1 € C[u,v], and s\,S2 to u- l,v ■ 1 respectively; 
since wp(2,3) = O(-b) and uj = Ox, 



degso = 5, and degsi, S2 = 7, 8. 



(4.7) 



The linear relations between so,si,S2 come from (fl.4j). Namely (essentially 
as in [ KM ], see also [PR], end of Section 1), sq,si,S2 correspond to the 



generators of L2 or the rows of M, and J M: L2 — > ii must map them to 
the image of F24 under the downarrow. This gives the equations 



or spelled out: 





(S2\ 








Si 








w 










.1 

y 



i 

.r 2 

-y 

—z 




z 

-yz 
-xz + y 1 


„2 



X 



t J 



I x 3 y + x 2 t \ 

4 2 
—x — xz 

xyz + zt 

-x 3 — z 2 

2x 2 y 

2 + ytJ 



V- 



-xy 



(4.8) 



I write R\, . . . , Rq for these relations. The equations f,g,h of T appear 
naturally 3x3 minors of J *M (or 2x2 minors of defective 2x3 blocks with 
zeros down one column) so that Ri, . . . , Rq can be solved by linear algebra 
to give expressions for (f,g,h) times S2,si,sq. For example 



fso 
hs 



o 00 

-3x yz + x y 



x 2 zt 



yt z 



x 6 + x 3 z 2 + xyzt + zt 2 , 



-x 5 y + 2x 2 yz 2 



xy 3 z + xz 2 t 



y 2 zt. 
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This makes explicit that so is a homomorphism 2r — > Ox- But it is much 
more useful to use R\, . . . ,Rq as they are, without taking determinants. 

5 The quadratic relations for sf , S\S2, s% 

Quadratic relations in sx,S2 must exist for several reasons. One is that the 
ring extension C[X] C C[Y] can be obtained by adjoining so first, followed 
by the normalisation C[X][so] C C[Y}; the second step is integral, and one 
guesses that C[Y] is generated by l,s\,S2 as a module over CLY][so]- Or, 
since s\,S2 correspond in some sense to uso and vs% on taking Poincare 
residue to T, and since, restricted to T, x = u 2 ,y = uv,z = v 2 + xu, there 
must be relations So, Si, S2 saying that 



S1S2 



xs 

■v4 



ZSq + XSqS\ 



G (xys ,tso,zs 2 ,x 2 z,xy 2 ,yt) , 
G (x2:so,y 2 so,a; 2 si,xyz,y 3 ,2:t} , 
€ (yzs ,xysi,tsi,x 4 ,xz 2 ,y 2 2:) , 



(5.1) 



where the right-hand side just lists alll monomials of degree 14, 15 and 16. 

To find these relations explicitly, I proceed as follows. Since MJ l M = 
and the first two rows of MN are zero by (|4.2|), the first two rows of M 



provide two syzygies between the 6 rows R\, . . . , Rq of (4.8). I massage the 
second very slightly to make them both 4-term syzygies: 



xR\ + yi?2 + zRz + tR^ = 0, and 
yRi + z(R 2 - xRi) + tR 5 + x 2 R G = 0. 



(5.2) 



This suggests realising the relations as 4 x 4 Pfaffians of 5 x 5 skew matrixess, 
as follows: 



Ix y 



\ 



z 

s 2 + x 2 



ys - xz 

-XS\ + ZSo 
X 3 ) 



so + y 

V 



X 
Si 
S 2 + X 1 



\ 



t 

■52 

-xy 

XS\ — ZSq f 



(5.3) 



My convention is to write only the upper triangular terms ai 2 , ai3, • • • , 045 of 
a skew matrix. The Pfaffians are defined by PUj.ki = a>ija>kl ~ a ika>ji + a u&jk 
(up to ±). 

The second matrix has first four Pfaffians R§, Rq, R2 — XR4, R±, and the 
fifth is one of the required quadratic relations: 



So 



2 2 1 

s 2 ~ zs + XSqSi 



yzs - 2xysi 



x 2 s 2 . 



(5.4) 
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The first matrix gives Ri, R 2 , R3, R4 and 

Pf 23.45 = xs\ - zsqSi + ys s 2 - x 2 ys - xzs 2 , 



(5.5) 



which is not quite what we want. However, subtracting S0-R5 gives a relation 
that is identically divisible by x, and taking out the factor gives 



So : s\ — xsl = xyso + zs 2 . 



(5.6) 



The equation for sis 2 can be found by constructing an explicit syzygy, 
for example by taking S1R4 + ySo + zR^, that has xs\s 2 as one term, and all 
terms of which are divisible by x. However, this is again a 4-term syzygy, 
once again best understood by realising the relations as Pfaffians: 



(x y 

z 

\ 



z 

Si 

s 2 + x 2 



si \ 

s 2 

—xso — 2xy 

si - yso J 



(5.7) 



Thus 



Si ■■ sis 2 - ysl = (y 2 - xz)so - x 2 si - 2xyz. 



(5.8) 



In fact, all the relations Ri, . . . , R5, So, Si, S 2 (but, annoyingly, not Re) can 
be written together as the 4x4 Pfaffians of the following 6x6 matrix: 



fx y 

z 



z 

Sl 



s 2 + x 2 



Sl 

s 2 

-xsq — 2xy 

-So - y s o 



\ 



yso - xz 
—xsi + ZSo 

x 3 

sqs 2 + x 2 s + x 2 y + xtj 



(5.9) 



6 Another example 

Consider the following two families of K3 surfaces: 

(a) The codimension 2 complete intersection C P(2, 2, 2, 3, 3), known 
to its friends as Fletcher 2 (14), with 9 x |(1,1) singularities on the 
plane P(2, 2,2) = P 2 . 

(b) The codimension 5 K3 surface Y C P(2, 2, 2, 2, 3, 3, 3, 3) in symmetric 
determinantal format, with 10 nodes. Write x±, . . . , X4, y±, . . . , for 
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coordinates, and let M be a general 4x4 symmetric matrix of linear 
forms in x%, . . . , x±. The ideal of Y is generated by the 14 equations 



M 



2/2 
2/3 

W 



0, y lVj = (Aft) 



here Aft is the adjoint matrix of signed maximal minors. Y is the 
determinantal hypersurface Y : (detM = 0) C P 3 = P(2,2,2,2), and 
the symmetric matrix M is the resolution of an ample divisorial sheaf 
Oy{A) (singular at the 10 nodes, where rankAf = 2), such that the 
Xi e H (O Y (2A), Vi G H (O Y (3A)). 

Then Y has a projection of Type IV from any node P to a codimension 2 
complete intersection Xq^ C P(2, 2, 2, 3, 3). 

Working top down from Y, I choose coordinates so that the node P = 
(1,0,0,0) in P(2, 2, 2,2), and 2/1,2/2 are local coordinates at P (compare 
[pPR|| , 3.4 for local coordinates on quotient singularities and their fractional 
divisor of zeros on a blowup). The geometric meaning of projection is that 
I blow up P to give an exceptional — 2-curve E, and consider the blown up 
surface polarised by the new Weil divisor A—^E. This eliminates x\ because 
it does not vanish on E\ moreover, it eliminates 2/1,2/2 because these have 
weight 3 in the graded ring, but vanish on E with multiplicity exactly i (see 
[pPR ], 3.4). This is the mark of a Type IV projection. It is an exercise to 



make the appropriate preliminary deduction about Af at P, and to do the 
elimination to get two equations of degree 6 in the remaining variables. 
Now work bottom up from X. As in Section |2[ one checks that 

P 1 -> P(2,2,2,3,3) given by (u, v) i-> (n 2 , uv, v 2 , u 3 , v 3 ) (6.1) 

(or a more general pair of cubics if you prefer) is an embedding, that the 
image is contained in the obvious conic q = x\x% — x\ and in the 3 sextics 
/ = x\ — 2/f, g = x\X2X^ — 2/12/2, h = x 3 — 2/2, and that two general sextic 
linear combinations of q, f, g, h define a quasismooth K3 surface ^6,6- 
I can resolve C[u, v] as a module over C\x, y] as in ([3.6]), obtaining 

C[«, v]< — L ^- Li J?- L 2 L 3 «- 0, (6.2) 
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where 



L = A®2A{-2), 
Li = 4A(-3) © 5A(-4), 
L 2 = 5A(-6) ©4A(-7), 
L 3 = 2A(-9) © A(-10). 

Here M is a skew 9x9 matrix with generic rank 6. Next, following the 
argument of (4.1-Q|), I can compose two maps between complexes to map 



the Koszul complex of X to 

A <- A(-6)ffiA(-6) <- A(-12) <- 

n i i (6-3) 

Lo <— L\ <— L2 <— <— 0, 

where the final downarrow defines an element col £ £2(12) analogous to 
Q4.5|) . The 3 unprojection variables X4, 7/3, 2/4 in degrees 2, 3, 3 correspond to 
the basis of L3, and as in ( |4.8p , the linear relations between them are in the 
Kustin-Miller form *P ^ p j = col. I believe that the quadratic relations for 

?/3, 2/32/4, 2/4 can be worked out as in Section |5[ 

This example is instructive: although both families of K3 surfaces are 
well understood, the equations describing the birational map between them 
are very subtle, and would require a lot of computation to elucidate. 



7 Magma routines 

7.1 Embedding Y and quasismoothness of X 4 

The first thing is to make the polynomial ring S = Q[u, v,x,y, z, t] and 
define the ideal of the graph of the embedding P(2, 3) — > P(4, 5, 6, 9): 

> Q := RationalsO; 

> S<u,v ) x ) y ) z,t> := PolynomialRing(Q, [2,3,4,5,6,9] ) ; 

> I := Ideal ( [-x+u~2,-y+u*v,-z+v~2+u~3,-t+u~3*v+v~3] ) ; 

(you can print / by doing "I;" at the prompt, and similarly throughout). 
Now eliminate the first two generators to get the ideal of the image T: 

> J := Eliminationldeal (I , 2) ; 

> [WeightedDegree(Basis(J) [i] ) : i in [1 . .#Basis(J)]] ; 
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answer: [20,18,19,21,22,23]. This says that the 1st, 2nd and 3rd basis 
elements have weight 20, 18, 19. (At this point I used the output of "J;" to 
write (2.2).) I make their generic linear combination of weight 24 by hand: 



> h := Basis(J)[l]; g := Basis(J) [3] ; f := Basis(J) [2] ; 

> F0 := x*h + y*g + z*f ; IsHomogeneous (F0) ; 

answer: true. 

Remark 7.1 Taking generic values such as 1, 1, 1 for coefficients in a linear 
system and jiggling them if necessary is a computer algebra substitute for 
Bertini's theorem. Quasismooth is an open condition, so if it ever holds, 
we would be infinitely unlucky to happen on a singular guy. In fact both 
F = xh + yg — zf and xh + yg give singular X24. 

To work in the polynomial subring R = k[x, y, z, t], the type checking in 
Magma insists that I force the elements x,y,z,t into it explicitly (otherwise 
it has no way of knowing that the user wants to identify the variables in 
R and S just because they have the same human names). One way is to 
set up a homomorphism (p: S — > R taking the generators u,v 1— ► 0,0 and 
x, y, z, t 1— > X, Y, Z, T: 

> R<X,Y,Z,T> := PolynomialRing(Q , [4,5 ,6 ,9] ) ; 

> fie := hom< S -> R I 0,0,X,Y,Z,T >; 

> F := fie(FO) ; 

Now I define the hypersurface X24 : (F = 0) C P(4, 5, 6, 9) and ask the 
big question: 

> PP:=Proj(R); X : =Scheme (PP , F) ; 

> IsNonSingular(X) ; 

answer: true. 

For completeness, I prove that P(2, 3) — > T is an isomorphism. For this, 
I check that monomials of weight i in C[X, Y, Z, T] map to polynomials in 
u,v that span the vector space of polynomials of weight i: 

> RR<U,V> := PolynomialRing(Q, [2,3] ) ; 

> psi := hom< R -> RR | U~2, U*V, U"3+V"2, IT3*V+V~3 >; 

> for i in [12. .24] do 

> K := ideaK RR | psi (MonomialsOf WeightedDegree (R, i) ) >; 

> [m in K : m in MonomialsOfWeightedDegree(RR,i)] ; 

> end for; 
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the last line returns the list of Boolean value (is ip(m) G K?) for every 
monomial in X, Y, Z, T of degree i G [12, ... ,24]. The answer [true, false, 
false], [true, true], etc., says that the 2 monomials U 3 and V 2 in degree 12 
fail the test, and every other monomial passes. 

7.2 Resolution of C[u, v] 

The problem is to treat C[u, v] as a module over the ring C[x, y, z, t], where 
the action is as in Q3.1| ). The following is based on a Magma routine written 
by Alan Steel and Gavin Brown. We treat C[u, v] as the quotient 

S/(x -u 2 ,y - uv, z-u? -v 2 ,t- u 3 v - v 3 ), (7.1) 

where S = C[u,v,x,y,z,t] (that is, the coordinate ring of the graph of i). 
It is certainly generated by l,u,v over the big ring S (in fact, by 1 only), 
and it is trivial to calculate the matrix of relations: 

/ u 

V 

o 

X — U 

y — uv 

S 2 
Z — U — V 

\t — u 3 v — v 3 

Define L C F to be the submodule of the free module F = S 3 generated 
by the rows of this matrix, and Grobnerise it. In Magma Grobner is a 
transitive verb that takes as its object a submodule of a free module. 

> Q := RationalsO; // Omit this line at your peril. 

> S<u,v,x,y,z,t>:=PolynomialRing(Q, [2,3,4,5,6,9] ) ; 

> Free := Module (S, 3) ; 

> L:=sub< Free | [ [u, -1, 0], [v, 0, -1], [-u~2+x, 0, 0], 

> [-u*v+y, 0, 0], [-u"3-v"2+z, 0, 0], 

> [-u"3*v-v"3+t, 0, 0] ] >; 

> Groebner (L) ; 

Grobnerising does many row operations on the matrix in ([7.2j), to put the 
dependence on u at the top and on the right, then the same for v. It gives 



-1 


o\ 





-1 























0/ 



0. 



(7.2) 
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the new basis of L: 



Basis (L) 



—y 





u\ 


—x 


u 





u 


-1 





— z 


X 


V 


-y 


V 





V 





-1 


xz + y 2 


x 2 





x 2 


—z 


y 


-z 2 


xz 


t 





xt 2 — y 2 z 2 





xt 


-yz 








-y 


X 


-zH 2 


y 2 z 3 


t 3 


yz 


-t 





-t 





z ) 



(7.3) 



V 



The rows not involving u, v generate the submodule N C F 
of all linear relations not involving u, v in the coefficients. 



S 3 consisting 



> LI := [f : f in Basis(L) | &and [Degree (f [i] ,u) le and 

> Degree (f [i] ,v) le : i in [1..3] ]]; 

> N := sub< Free I LI >; 

> MinimalBasis(N) ; 



answer: the minimal basis of ./V is the matrix (|3.3j ), up to cosmetic changes 



7.3 The K3 Database in Magma 



Export 2.8 of the computer algebra system Magma [Ma] includes the first 
public version of Gavin Brown's K3 surface database, which is largely a 
reworking in computer language of the calculations of Anthony Iano-Fletcher 
and Selma Altmok's Warwick theses |Q, Q. It lists families of polarised 
K3 surfaces in searchable form. Further mathematical presentation, and 
more computer documentation will be available soon (see [ABR]). Here I 
only illustrate how to use the database to find the projections of Type IV 
mentioned in the introduction. 

> D:=K3Database("t") ; 

> #D; 
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answer: the database D contains 391 elements. The first line defines D to be 
the K3Database, and requests that Magma refers to the variable in Hilbert 
series by its human name t, rather than its preferred internal representation 
$.1. (You can put empty brackets () if you like, but if you omit the brackets, 
Magma will just sulk and refuse to play.) 

> #[ X : X in D | Codimension(X) eq 4 ] ; 

answer: 142. Note the logic of Magma's very useful and flexible sequence 
construction [ ] . Here it constructs the sequence of X € D such that 
codim V = 4. It runs through all X G D (the entry after the colon), evalu- 
ates the Boolean value of the statement following the bar | (is codimX = 4?), 
and returns X (the entry before the colon) if the answer is yes. We could 
ask it to return any function of X by putting f(X) before the colon. Asking 
for the number of items in a list # [ ] is usually preferable to asking for the 
list itself, which may be huge and uninformative. An alternative is to give 
the list a name, and ask to see just a bit of it: 

> cod4 := [X : X in D | Codimension(X) eq 4 ] ; 

> cod4[l. .5] ; 

To play with projections, we first ask Magma to calculate the possible 
projection of Type I, II and IV of each element of the database D from each 
of its singularities: 

> Centres (~D) ; 

> Cod4Type4 := [X : X in D | Codimension(X) eq 4 \ 

> and &or[ p[2] eq [4] : p in Centres (X)] ]; 

The last line here defines Cod4Type4 to be the sequence of X £ D having 
codim X = 4 and such that at least one of the centres p of X has Type IV. 
The internal [ ] constructs a sequence of Boolean values (for p a singularity 
of X, is the projection from p of Type IV?) and the &or[ ] evaluates "or". 
To write this code you need to know a little about how the database stores 
its data. You can find out by interrogating it. For example, 

> X := D[293]; Centres (X); 

tells you that No. 293 in D has a singularity of type |(2, 3) with a projection 
of Type I, etc., and you figure out that each centre p is stored as a sequence of 
entries of which the second p[2] is the type of the projection. I ask how many 
elements does Cod4Type4 contain, and what are the first two elements? 
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> #Cod4Type4; 

> Cod4Type4[l. .2] ; 



answer: 5 elements. For the first, see below. The second is the example 
treated in Sections [2]-|5] above. All five of these K3s are from Selma Altmok's 
thesis. For historical reasons, the database lists the codimension 4 K3s as 
Altinok4(-24) through to Altinok4(121), missing Nos. 0, 76, 81, 107. The 
5 elements of Cod4Type4 are listed as Altinok4(84), (111), (53), (-11) and 
(-6). 

Example 7.2 The first element of the above list is Altinok4(84), the codi- 
mension 4 K3 surface Y C P(5, 5, 6, 7, 8, 9, 11) with basket 2 x |(2, 3), ^-(5, 6) 
and Hilbert numerator 

1 _ t u - t 15 - 2t 12 - t 17 - 2t 18 -t 19 -t 20 + --- 

The database says that this has a projection of Type IV to a surface hav- 
ing the same numerical type as Xqq C P(5, 6, 22, 33), but containing a 
copy of P(2,3). This could happen because Xqq is actually a hypersur- 
face of P(5, 6, 22, 33) containing the image of a general embedding P(2, 3) «— > 
P(5, 6, 22, 33); or Xqq could itself be a bit degenerate, for example, with the 
curve P(2, 3) as a base locus of some monogonal linear system. It would be 
interesting to see this and the other cases in the list worked out in the spirit 
of this paper. 



Remark 7.3 We expect subsequent exports of Magma to contain substan- 
tial improvements to the internal structures of the K3 database, its contents, 
and its user interface. Whereas at present we use Magma mainly as a book- 
keeping device to handle large quantities of combinatorial data, we hope 
that in future we can use it to automate many of the computations in com- 
mutative and homological algebra involved in the geometry of K3 surfaces 
and Fano 3-folds, including the calculations above and those in [CPR] and 



[Ki 
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